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It is shown that dispersion of waves in the magnetosphere regarded as a cold plasma may con- 
veniently be handled from whistler frequencies down to frequencies small compared with the ionic 
gyrofrequency by an approximation (22) similar to the Appleton-Hartree approximation (23) used at 
radio frequencies. In both cases the ratio of the ionic mass to the electronic mass is assumed infinite, 
but whereas for the radio approximation this is achieved by taking the ionic mass infinite, it is achieved 
for the hydromagnetic approximation by taking the electronic mass zero. This permits extension of 
the concept of quasi-longitudinal and quasi-transverse approximations to the entire frequency spec- 
trum (fig. 2); propagation is predominantly quasi-transverse at frequencies small compared with the 
ionic gyrofrequency. The directional pulse patterns for radiation of a pulse from a source are cal- 
culated (fig. 11) and the transition from the omnidirectional MHD wave to the whistler mode is illus- 
trated (fig. 13). The hydromagnetic approximation here used facilitates detailed calculation of 
elliptical polarization (figs. 14-19) and space-charge (fig. 20). Attenuation due to ionic collisions with 
neutral gas is calculated (fig. 21) on the assumption that the neutral gas is too massive to take part 
in the wave motion. 



1. Introduction 

For many purposes the ionosphere and the magnetosphere may be regarded as a cold mag- 
netoplasma. Considerable interest attaches to understanding the propagation of waves in such 
a plasma from frequencies small compared with the ionic gyrofrequency to frequencies large 
compared with the electronic plasma frequency. The radio frequency end of the spectrum has 
been studied for a number of decades and is understood in considerable detail. Reference may 
be made to books by RatclifTe [1959] and by Budden [1961]. The hydromagnetic end of the spec- 
trum was developed independently by Alfven [1950], and has become an important aspect of fluid 
mechanics. It has long been realized, however, that, if both ions and electrons are taken into 
account, the phenomena at the two ends of the spectrum are both incorporated in the magneto- 
ionic theory [Astrom, 1950]. The unity of the theory at radio and hydromagnetic frequencies has 
been brought out in recent books by Denisse and Delcroix [1961], by Stix [1962], and by Allis, 
Buchsbaum, and Bers [1963]. However, the complete version of the magneto-ionic theory required 
to cover the entire spectrum is algebraically complicated. For radiation from a source in a cold 
uniform magnetoplasma, the various topological forms taken by the phase propagation surfaces 
have been ascertained and are conveniently demonstrated in a diagram of the type introduced by 
Clemmow and Mullaly [1955], and modified by Allis [1959]. However, by making appropriate 
approximations, one can concentrate on those situations of greatest interest in the ionosphere and 
the magnetosphere, and this permits a simpler treatment that facilitates greater understanding. 
The best known simplification of this type is, of course, the Appleton-Hartree formula, [Appleton, 
1932; Hartree, 1931] together with its quasi-longitudinal and quasi-transverse approximations 
[Booker, 1935]. We shall describe the Appleton-Hartree formula as based on the radio approxi- 
mation to the complete magneto-ionic theory. It is the objective of this paper to make cor- 
responding use of a similar approximation covering the hydromagnetic end of the spectrum. 
Together, the hydromagnetic and radio approximations cover the vast majority of situations of 
practical interest in the magnetosphere regarded as a cold plasma. 

In both the radio and the hydromagnetic approximations it is assumed that the ratio of the 
mass of an electron to the mass of an ion is negligibly small. In the radio approximation this is 
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achieved by allowing fully for the inertia of the electrons but taking the mass of an ion as infinite. 
In the hydromagnetic approximation, on the other hand, we allow fully for the inertia of the ions, 
but take the mass of an electron as zero. At whistler frequencies, under conditions of applicabil- 
ity of the Eckersley law [Eckersley, 1935], the two approximations give the same results, thus 
permitting a good fit between the hydromagnetic and radio approximations under conditions of 
practical interest. 

The two elliptically polarized waves under discussion are known at radio frequencies as the 
ordinary and extraordinary waves (O and X waves). The O wave is the wave for which propagation 
becomes independent of the imposed magnetic field when the direction of this field is turned so 
as to be perpendicular to the direction of phase propagation. 3 Below the ionic gyrofrequency 
these waves may also be referred to as the slow and fast waves respectively. At frequencies small 
compared with the ionic gyrofrequency they are referred to in the fluid mechanics literature 
[Bershader, 1959] as the transverse and fast waves, and the term slow is reserved for the acoustic 
mode. The two waves may also be classified according to the direction of rotation of the electric 
vector round the direction of the imposed magnetic field. This is actually the most convenient 
classification within the confines of the hydromagnetic approximation used in this paper. In this 
approximation the assumption of zero electronic mass gives infinite conductivity in the direction 
of the imposed magnetic field. This makes the plane of the electric ellipse perpendicular to the 
imposed magnetic field for all directions of propagation. We shall describe waves as left-handed 
or right-handed (L or R waves), according to whether the elliptical polarization of the electric 
vector is left-handed or right-handed about the direction of the imposed magnetic field. We 
shall draw diagrams for the L and R waves, respectively, but particularly for the R waves, we 
shall indicate which parts of the diagram refer to the O wave and which to the X wave. 



2. Dispersion Formulas 

We shall consider wave propagation in a cold magnetoplasma consisting of electrons of mass 
m e and charge q e vibrating with velocity v e , and ions of mass mi and charge q\ vibrating with velocity 
\i. We shall neglect nonlinear effects and electronic collisions. Initially we shall neglect ionic 
collisions, but in section 7 allowance will be made for ionic collisions with neutral particles in 
circumstances when the neutral gas is too massive to take significant part in the wave motion. 
Let B be the imposed magnetic flux density and E the electric field of the wave. In the absence 
of collisions the equations of motion for the electrons and ions are 

m e -^ = qe(E + v e XBo) (1) 

Tm^=qi(E+ViXBo) (2) 

where t denotes time. The radio approximation is based on the assumption that nn is infinite, 
so that ionic motion is unimportant and (2) disappears; on the other hand, no approximation is 
made in (1). The hydromagnetic approximation used in this paper is based on the assumption 
that m e vanishes, so that (1) reduces to 

= E+v,xB . (3) 

On the other hand no approximation is made in (2). Equation (3) implies that an electron moves 
in such a way as to keep the total force on it zero. For a finite electronic velocity parallel to the 
imposed field there can be no component of E in this direction. The electric field of the wave 
therefore vibrates in planes perpendicular to the imposed magnetic field. 



3 Note that Astrom [1950] reverses the names "ordinary" and "extraordinary" at hydromagnetic frequencies. 
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To describe the propagation of a wave in the plasma the equations of motion of the electrons 
and ions must be combined with Maxwell's equations and with the equations of continuity. Except 
for propagation strictly parallel to the imposed magnetic field, waves even in a cold magnetoplasma 
involve compression and rarifaction of both the electrons and the ions. 4 If N e and Nt are the mean 
electron density and ion density respectively, and 8N e and 8/V, are the corresponding fluctuations 
caused by the wave, the equations of continuity are 

ddN e M 

ir—wv-*. ( 5) 

The current density J and the charge density p in the plasma are given by 

J = N e qe\ e + NiqiYi (6) 

and 

p = q e 8N e + qi8Ni. (7) 

It is assumed that the plasma is, on the average, neutral so that 

0=q e N e + qiNi. (8) 

The equation of continuity of charge need not be written down since it is derivable from (4) and 
(5) by use of (6) and (7). 

To obtain the dispersion equation for a wave of angular frequency co and propagation vector k 
we substitute into the above mentioned equations the wave function 

exp {j(<Dt — k- r)} (9) 

where r denotes position vector. The equations of motion (3) and (2) then become 

= E + v,xB o (10 ) 

joomiVi = qt(E + vi X B ) (11) 

while the equations of continuity become 

M e 1 

~/v7 = * k,v * (12) 

— -^- (13) 

On eliminating the magnetic field of the wave from Maxwell's equations we obtain 

(k • E)k- ( k 2 - ^j- J E=JG)fJL (N e qeVe + NiqiYi). (14) 



4 Thus equation (101) of Lighthill's treatment [1960] is inapplicable to a cold magnetoplasma and this upsets the corresponding section of his paper except for 
propagation strictly parallel to the imposed magnetic field. 
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When required, the magnetic field H of the wave is given by 

H=— kxE (15) 

where fio is the vacuum inductivity. Equation (15) shows not only that the electric and magnetic 
fields of the wave are mutually perpendicular but also that the magnetic field of the wave vibrates 
in a plane perpendicular to the direction of phase propagation. 

In the magnetosphere we are interested in discussing the propagation of waves above and 
below the following frequencies: the electronic plasma frequency /v e , the electronic gyrofrequency 
JMe, the ionic plasma frequency /m, and the ionic gyrofrequency /mi- The corresponding angular 
frequencies will be denoted by com?, u>Me, com and com?. At hydromagnetic frequencies it is conven- 
ient to express velocities in terms of the ionic Alfven velocity A\ given by 

4i = c— (16) 

(*>Ni 

where c is the velocity in free space. At whistler frequencies it is convenient to compare velocities 
with the electronic Alfven velocity A e given by 

A e = c—- (17) 

Let the imposed magnetic field point in the direction of the z-axis, and let the direction of phase 
propagation lie in the yz plane as shown in figure 1. With this choice of coordinate axes we may 
write 

Bo = (0,0,B ) (18) 

k = (0, k ± , k\\) = k{0, sin 6 P , cos 6 P ). (19) 

The angle 6 P is the angle between the direction of phase propagation and the direction of the im- 
posed magnetic field. It must be distinguished from the angle 6 g introduced later for the angle 
between the direction of group propagation and the direction of the imposed magnetic field. We 



X AXIS DOWN 

FIGURE 1. Illustrating the direction of the imposed mag- 
netic flux density B and the direction of the phase 
propagation vector k in relation to the axes of 
coordinates. 
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shall encounter situations in which the direction of phase propagation is perpendicular to the 
imposed magnetic field while the direction of group propagation is parallel to the imposed magnetic 
field. 

If we substitute from (18) and (19) into (10) through (14), we obtain a set of linear algebraic 
equations for the ratios of the following quantities: 



8N e > 8Ni 9 E x , E y , Vex, V ey , Vez, V ix , V iy . 



(20) 



The component of the ionic velocity in the direction of the imposed magnetic field does not appear 
because there is no component of the electric field in this direction to drive the ions; thus the 
ionic motion, which is substantially the fluid motion, is in planes perpendicular to the imposed 
magnetic field regardless of the direction of propagation. Elimination of the quantities (20) from 
the linear algebraic equations gives the dispersion relation 
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of which the solution is 
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(22) 



This is the formula that, at hydromagnetic frequencies, replaces the Appleton-Hartree formula 



kc 




(23) 



i 2 a mK 

, "+COS z 6 P r 



applicable at radio frequencies. The alternative signs in (22) and (23) refer to the two magneto- 
ionic waves. 

The quantity kc/co appearing in (21), (22), and (23) is the refractive index referred to the free 
space velocity c in the manner customary at radio and optical frequencies. However, at hydro- 
magnetic frequencies, it is often more convenient to use the refractive index kAi/co referred to the 
ionic Alfven velocity. Furthermore, throughout the magnetosphere and in many other situations 
encountered in plasma physics, the ionic gyrofrequency is small compared with the ionic plasma 
frequency. Consequently the expression 



2 

i I Mm 

9 



cot 



(24) 



appearing in (21) and (22) may usually be replaced by unity. Hence (21) and (22) may usually be 
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replaced by 



and 



^('-^)(^) , -«™.>('-|)(^)* + ('-^°)-<— » 
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cos 2 # p [2 * \ w^, 
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^VI^FIF^I}] 



+ cos 2 P ^- [ | (o)Mi < com). (26) 



Moreover, at frequencies of the order of the ionic gyrofrequency and below, (26) may be further 
simplified to 



ffl-t 



1 + — rz"lo sin2 °p 

cos z Op [Z 



"m 



1 CO 2 
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(co and comi <^ com). (27) 



When the frequency becomes small compared with the ionic gyrofrequency, the dispersion formula 
(27) tends to the one encountered in fluid mechanics literature on magnetohydrodynamic waves. 
Virtually all the techniques that have been developed for handling the radio approximation 
(23) are equally applicable for studying the hydromagnetic approximation (22). In particular, 
the concept of approximations for quasi-longitudinal and quasi-transverse propagation that have 
proved convenient in handling (23) is equally applicable to (22). We have quasi-longitudinal 
propagation when the cos 2 P term under the radical dominates the sin 4 P term; we have quasi- 
transverse propagation when the sin 4 P term dominates the cos 2 6 P term. By equating these two 
terms in (22) we see that, for the hydromagnetic approximation, the transition from quasi- 
longitudinal to quasi-transverse propagation occurs where 



sin 2 Op 
2 cos Op 



CO 

CO\ti 



1 + 



<»Ni 



CO 
CO? 



(28) 



If expression (24) can be replaced by unity, (28) becomes 



sin 2 Op 


CO 
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<4i 



(comi <^ com). 



(29) 



Other useful approximations to (28) are 



sin 2 ft 



'£-= ^Ni 



and 



2 COS Op COCJMi 



sin 2 D 



2 COS Op COMi 



(co p> com and com) 

(co and co M i <com). 
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(30) 



(31) 



Figure 2. Illustrating the distribution between quasi- 
longitudinal and quasi-transverse propagation in the 
magnetosphere roughly one earth's radius out in the 
equatorial plane. 
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For the radio approximation (23) the transition from quasi-longitudinal to quasi-transverse prop- 
agation occurs where 



sin 2 O p 
2 cos dp 



1 <»Ne 
0> 2 



a> 



(32) 



to which the low frequency approximation is 



sin' 



-&L-. 



4e 



2 COS Op (OCOMe 



((D <^ (t) Ne ). 



(33) 



It should be noted that the expressions quasi-longitudinal and quasi-transverse refer to the direc- 
tion of propagation and not to the direction of vibration of particles. It should further be noted 
that these expressions refer to the direction of phase propagation in relation to the imposed mag- 
netic field and not to the direction of group propagation. Situations will occur in which propa- 
gation is quasi-transverse even though the direction of group propagation is nearly parallel to the 
imposed magnetic field. 

To obtain some idea when propagation is quasi-longitudinal and when quasi-transverse, let 
us apply (28) through (33) to a magnetoplasma of the type existing in the magnetosphere roughly 
one earth's radius out in the equatorial plane. The frequencies /mu /nu fiue, and f^ e then have 
values such as those shown in figure 2. The continuous curve in figure 2 is drawn from (28) and 
the broken curve from (32). We notice how well the hydromagnetic approximation joins onto the 
radio approximation. The reason is that (30) and (33) are numerically identical because the ratio 
of the square of the plasma frequency to the gyrofrequency is independent of particle mass. 
Figure 2 illustrates the well-known fact that, in the radio range, propagation is mainly quasi-longi- 
tudinal except near the electronic plasma frequency. Figure 2 also shows that, in the hydro- 
magnetic range, propagation is principally quasi-transverse below the ionic gyrofrequency and 
principally quasi-longitudinal above. The usual dispersion equation for magnetohydrodynamic 
waves given in the fluid mechanics literature is in fact a quasi-transverse approximation, and 
this is why such papers sometimes make no reference to circularly polarized waves. 

We have seen in figure 2 that the hydromagnetic approximation fits on to the radio approxi- 
mation remarkably well, and we need to investigate to what extent this convenient result is generally 
applicable. For frequencies high compared with the ionic gyro and plasma frequencies (22) 
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may be written 
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while, for frequencies small compared with the electronic plasma frequency, (23) may be written 
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(34) and (35) would be completely identical were it not for the term (i> 2 yJ Ja) 2 Me in the denominator in 
(35). Unfortunately, in the magnetosphere, the electronic plasma frequency is not small compared 
with the electronic gyrofrequency, and so the additional term in (35) as compared with (34) is 
significant. However, this is only true for quasi-transverse propagation, which, as indicated in 
figure 2, is comparatively unimportant at the frequency where the transition is made from the 
hydromagnetic approximation to the radio approximation. For quasi-longitudinal conditions at 
the transition frequency the cos 2 d p term under the radical in (35) dominates the sin 4 d p term, so 
that (34) and (35) are identical if the transition frequency is small compared with the electronic 
gyrofrequency. In fact, for quasi-longitudinal propagation, (34) and (35) reduce to 



-Y = i± — S&-— (QL) 

to / (DtoMi COS Op (37) 



kcV 

(O ) <t)<i>Me COS dp 



= 1± ,..,JT^ (QD (38) 



and the identity of these equations follows from (36). It is upon (38) that the Eckersley's law for 
whistlers is based. Equations (37) and (38) show that, while the basis of the Eckersley law is 
usually thought of as the low frequency behavior of the radio approximation, it may equally well 
be thought of as the high frequency behavior of the hydromagnetic approximation. Thus in all 
circumstances where there is an interval of frequency in which the Eckersley law for whistlers is 
applicable, the hydromagnetic approximation is bound to fit smoothly onto the radio approximation. 
The fit between the hydromagnetic and radio approximations is illustrated in figure 3 for 
longitudinal propagation of the R wave and for the same ionization density and imposed magnetic 
field used in figure 2. The continuous curve in figure 3 shows the dispersion behavior for the R 
wave using the full cold plasma theory. The broken curve in figure 3 shows the Eckersley law, 
and we see how this agrees with the full theory at audio frequencies. If the dispersion curve were 
plotted entirely on the basis of the radio approximation it would follow the continuous curve in 
figure 3 from whistler frequencies upwards but the broken curve from whistler frequencies down- 
wards. On the other hand if the dispersion curve were plotted entirely from the hydromagnetic 
approximation it would follow the continuous curve in figure 3 from whistler frequencies down- 
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Figure 3. Illustrating the fit between the hydro magnetic 
and radio approximation for longitudinal propagation 
of the R wave about one earth's radius out in the equa- 
torial plane. 
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FIGURE 4. Illustrating propagation of the X wave trans- 
verse to the earth's magnetic field about one earth's 
radius out in the equatorial plane. 



wards but the broken curve from whistler frequencies upwards. Thus figure 3 illustrates for longi- 
tudinal propagation how the hydromagnetic approximation fits onto the radio approximation at 
whistler frequencies where both approximations give the Eckersley law. 

On the other hand, when the direction of phase propagation is transverse to the imposed 
magnetic field, there is misfit between the hydromagnetic and radio approximations at whistler 
frequencies on account of the additional term in the denominator of (35) not present in the denomi- 
nator of (34). This misfit is illustrated in figure 4 for the ionization density and imposed magnetic 
field used in figure 2. Figure 4 shows that, for propagation of the X wave transverse to the imposed 
magnetic field, there is a stop-band extending from the lower critical frequency in the radio range 
down to a frequency of/*, close to the ionic plasma frequency. The misfit occurs in this stop-band 
and is indicated by the broken curve joining the two approximations. It should be mentioned 
that the lower edge of the stop-band in figure 4 does not occur exactly at the ionic plasma frequency. 
Furthermore, the effect of the misfit between the hydromagnetic and radio approximations shifts 
the low frequency limit of the stop-band downwards to some extent; indeed, far out in the magneto- 
sphere, the shift is as much as a decade. It is feasible to modify the hydromagnetic approximation 
to avoid this misfit, but this modification is not discussed in this paper. 

The hydromagnetic and radio approximations cover almost everything required for propa- 
gation in the magnetosphere regarded as a cold plasma, except for quasi-transverse propagation 
at whistler frequencies. This exception is quite restricted as shown in figure 2. Virtually every- 
thing of interest in the cold plasma theory of wave propagation in the magnetosphere is therefore 
obtained by combining the Appleton-Hartree formula at radio frequencies with the hydromagnetic 
approximation given by (22). 



3. Dispersion Curves 

To illustrate the frequency behavior of the hydromagnetic dispersion formula (22), it will be 
necessary to assume a numerical value for the ratio of ionic plasma frequency to ionic gyrofre- 
quency. In the magnetosphere this ratio varies from about 50 to 1,000. Use of too large a value 
of this ratio is inconvenient, however, as some of the phenomena to be illustrated occur for values 
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of 6 P between tan 1 0Jml(0Mi and \ir. We shall therefore use for purposes of illustration the value 

com c 



0)Mi 



A, 



= 10.7 



(39) 



which corresponds to 



(*>Ne _ 
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fV4 (hydrogen) 
Vs (helium) 
[Vi6 (atomic oxygen). 



(40) 



However, the curves will be plotted so as to be independent of the value of the ionic plasma fre- 
quency for wave frequencies of the order of the ionic gyrofrequency and below. This is achieved, 
in accordance with (27), by referring the refractive index not to the free space velocity but to the 
ionic Alfven velocity. 

For the numerical values given in (39) and (40), the transition between quasi-longitudinal and 
quasi-transverse propagation is illustrated in figure 5 on the basis of (28). In figure 5a eq. (28) 
may in fact be replaced by (31), and in figure 5b by (29). Figure 5b illustrates the wide range of 
conditions for which propagation is quasi-longitudinal. Moreover, as the ratio of the ionic plasma 
frequency to the ionic gyrofrequency increases above the value given in (39), the regions of quasi- 
transverse propagation above the ionic gyrofrequency become even more restricted than is shown 
in figure 5b. On the other hand, for frequencies of the order of the ionic gyrofrequency and below, 
the behavior is independent of the value of the ionic plasma frequency as shown by (31) and figure 
5a. As the frequency tends to zero, propagation becomes predominantly quasi-transverse no 
matter how large the value of (omIwmu 
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Figure 5. Illustrating for the hydro magnetic approxi- 
mation the transition from quasi-longitudinal to quasi- 
transverse propagation (a) on the scale of the ionic 
gyrofrequency and (b) on the scale of the ionic plasma 
frequency. 
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To study the variation of refractive index with frequency it is convenient initally to examine 
the situations in which the refractive index for one of the waves is either zero or infinite, or in 
which the refractive indices for the two waves are identical. The zeros of refractive index are 
obtained by equating to zero the term in (21) not involving kc/co. The frequency derived is outside 
the range of validity of the hydromagnetic approximation under conditions of interest in the mag- 
netosphere. Zeros of refractive index occur in the range of validity of the radio approximation 
and are given by the well-known ionospheric formulas for critical frequency. For the magneto- 
sphere there are no zeros of refractive index at hydromagnetic frequencies. 

Besides the infinity of refractive index occurring at radio frequencies, there are infinities 
occurring at hydromagnetic frequencies. Infinities of refractive index for the hydromagnetic 
approximation are obtained by equating to zero the coefficient of (kc/co) 4 in (21). We see that 
there is an infinity of refractive index when co = comx and when d p = 1 /27T. Thus, for the hydro- 
magnetic approximation, the refractive index of one of the waves is infinite for all frequencies when 
the direction of phase propagation is transverse to the imposed magnetic field. This is in fact 
the O wave for which transverse propagation is independent of the imposed magnetic field. For 
transverse propagation the O wave has a stop-band extending from the electronic plasma frequency 
indefinitely downwards. At hydromagnetic frequencies this wave is highly evanescent and has 
a refractive index on the imaginary axis close to the point at infinity. This becomes identical 
with the point at infinity for the hydromagnetic approximation because of the neglect of electron 
inertia. Thus, by putting P = V27T in (21), we see that the refractive indices of the O and X waves 
for strictly transverse propagation using the hydromagnetic approximation are 



and 
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(T, X). (42) 



It is from (42) that the hydromagnetic approximation shown in figure 4 is plotted. For trans- 
verse propagation of the X wave there is an infinity of refractive index for the hydromagnetic 
approximation at the frequency /«, (angular frequency co x ) for which the denominator in (42) van- 
ishes. Hence 

(0oo = (<O 2 m + <ti 2 m )V2 (43) 

which reduces to co x = com when comi ^ com- As mentioned in connection with figure 4 misfit be- 
tween the hydromagnetic and radio approximations for quasi-transverse propagation at whistler 
frequencies depresses the value of co x below that given by (43). When the effect of electron inertia 
at the ionic plasma frequency is taken into account Stix [1962] gives for co^ the approximate equation 

— = — i— + — — 

to 2 * to 2 Ni + to 2 m COMiCOMe ^ ' 



which may usually be written in the form 

1 1 



L)2 C0 2 m COMiCOMe 



{coMi < com). (45) 



For the curves drawn in this paper the frequency co x will be indistinguishable from com, but it should 
be remembered that co x is given more generally by (45). The difference is numerically important 
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at the higher levels in the magnetosphere and also in laboratory plasma physics. This feature has 
been studied by C. O. Hines [1957] and by Auer, Hurwitz, and Miller [1958]. 

Before actually drawing dispersion curves it is important to notice the conditions under which 
the refractive indices for the two waves are identical. For the hydromagnetic approximation this 
happens when the quantity under the radical in (22) vanishes, that is, when a) = w x and 6 p = 1 /2 it. 
This is simply the situation in which the expression on the right-hand side of (42) is infinite and is 
consequently identical with that on the right-hand side of (41). As is well known for the radio 
approximation, situations where the refractive indices of two waves are identical lead to coupling 
between the waves and to interchange in the names of the waves. This behavior for the hydro- 
magnetic approximation is conveniently illustrated by plotting the inverse of the square of the 
refractive index versus frequency as shown in figure 6. The continuous and broken curves for 
P = 9O° in figure 6 are based on (41) and (42), respectively and we notice that they intersect near 
the ionic plasma frequency. As we depart from strictly transverse propagation, the curves for 
the O wave become connected to the curves for the X wave in the manner illustrated. The sharp 
change as we move away from strictly transverse propagation conditions becomes a continuous 
but rapid change if electron inertia is taken into account or if ionic collisions are taken into account. 
The behavior is qualitatively similar to that encountered in the radio approximation close to strictly 
longitudinal propagation conditions at the electronic plasma frequency [Booker, 1934]. 

We are now in a position to replot figure 6 to give the dependence of refractive index upon 
frequency. This is done in figure 7. From figure 7b we see that the R wave is the X wave for 
/</«, and the wave for/>/ x . From figure 7a we see that the L wave has a hydromagnetic 
pass-band extending from zero frequency to the ionic gyrofrequency and that, in this band, the 
L wave is also the O wave. Near the low frequency end of this band the L wave is the wave 
originally introduced by Alfven [1950], while near the high frequency end of the pass-band the 
L wave is what is known to laboratory plasma physicists as the ion cyclotron wave [Stix, 1962]. 
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Figure 6. Illustrating the exchange in the nomenclature 
of the O and X waves near the ionic plasma frequency. 



Figure 7. Illustrating, for the hydromagnetic approxi- 
mation, refractive index versus frequency (a) for the L 
wave and (b) for the R wave. 
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Various sections of the dispersion curves shown in figure 7 (and of subsequent curves) may 
be conveniently approximated by using the quasi-longitudinal and quasi-transverse approximations. 
When the cos 2 6 P term under the radical dominates the sin 4 P term, (22) reduces to the quasi- 
longitudinal approximation 
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where the upper sign refers to the L wave and lower sign to the R wave. By putting 6^ = 0, we 
deduce for strictly longitudinal propagation using the hydromagnetic approximation 
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When the sin 4 P term under the radical in (22) dominates the cos 2 6 P term, we obtain for the quasi- 
transverse approximation for the O wave 
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of which useful further approximations are 
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To obtain the quasi-transverse approximation for the X wave the radical in (22) must be expanded 
by the binomial theorem, and we obtain 
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which may also be written as 
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From (48) and (52) we easily verify the truth of (41) and (42) when the direction of phase propaga- 
tion is strictly transverse to the imposed magnetic field. 
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As an example of the utility of the quasi-longitudinal and quasi-transverse approximations, 
consider the curve for 6 P = 60° in figure 7b. From figure 5b we see that, over most of the relevant 
frequency range, propagation is quasi-longitudinal. Thus much of the 60° curve in figure 7b is 
represented by (46) using the lower sign. On the other hand, for the 88° curve in figure 7b, we 
see from figure 5b that propagation passes from quasi-transverse at sufficiently low frequencies 
to quasi-transverse at sufficiently high frequencies through an intermediate frequency range of 
quasi-longitudinal propagation. The corresponding sections of the 88° curve in figure 7b may 
therefore be sketched in from (51), (48), and (46) using the lower sign. In particular, the value to 
which the 88° curve in figure 7b is leveling off as the frequency moves above the ionic plasma 
frequency is given by (50). 



4. Group Propagation 

We now wish to examine group propagation for the hydromagnetic approximation and to see 
how it fits in with group propagation for the radio approximation. We are here considering a 
situation in which a pulse whose spectrum is principally confined to a limited band of frequencies 
is radiated into a homogeneous magnetoplasma from a source of limited dimensions near the 
origin. The radiated pulse may be regarded as an aggregate of wave packets radiated from the 
origin in different directions with different group velocities. At any time subsequent to radiation 
these wave packets lie along a surface that constitutes the directional pulse pattern of the source 
in the magnetoplasma. We wish to understand how the directional pulse patterns for both mag- 
netoionic waves vary with frequency from radio frequencies down to frequencies small compared 
with the ionic gyrofrequency. 

Application of the principle of stationary phase to the wave function (9) gives 

r da) 

7 = ak (53) 

which shows that the wave packets travel radially outwards from the origin with the group velocity 

9 5k (54) 

We shall plot polar diagrams of the group velocity vector V^ for both magnetoionic waves and for 
a series of frequencies. Equation (54) shows that, if the dispersion equation is regarded as de- 
fining in k space surfaces of constant o>, then the gradient of the contour markings on these surfaces 
gives the group velocity in magnitude and direction. For this purpose it is convenient to write 
the dispersion equation (21) as 



V V W \ u>% ( <a%,J co 2 M .\(om/ IA 
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(55) 



If we regard (k, 6 P ) as polar coordinates, (55) gives a curve for each value of w. These curves 
are shown in figure 8 for the L wave and in figure 9a for the R wave. If we take the mirror image 
of these diagrams in the axes and rotate about the imposed magnetic field, we obtain surfaces 
with contour markings whose gradient gives the group velocity vector in accordance with (54). 
In figure 8 the transition between quasi-longitudinal and quasi-transverse propagation is shown, 
and we notice that the arrangement of straight parallel contours familiar in the fluid mechanics 
literature concerning magnetohydrodynamic waves is restricted to the quasi-transverse region. 
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FIGURE 8. Illustrating, for the L wave, surfaces of constant 
w in k space using the hydromagnetic approximation. 
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Figure 9. Illustrating, for the R wave, surfaces of constant co in k space (a) using the hydro- 
magnetic approximation on the scale of the ionic plasma frequency and (b) using the radio 
approximation on the scale of the electronic gyrofrequency. 



In figure 9a the broken contours, unlike the continuous contours, form oval curves when the other 
three quadrants are drawn. In figure 9a the propagation is quasi-longitudinal over the entire 
diagram except for a narrow strip along the perpendicular axis. The way in which the hydro- 
magnetic approximation for the R wave shown in figure 9a fits onto the radio approximation may 
be seen from figure 9b. Figure 9b illustrates the whistler dispersion formula 
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and is based on the work of R. L. Smith [1960] and R. Gendrin [I960]. Whereas figure 9a is drawn 
on the scale of the ionic plasma frequency, figure 9b is drawn on the scale of the electronic gyro- 
frequency. To fit the two diagrams together the contour map in figure 9a must be appropriately 
contracted in scale and fitted between the perpendicular axis and the 0.1 fuje contour in figure 9b. 
We use the contour maps in figures 8 and 9 in the following way. If we are interested in 
radiation of a wave packet from the origin in a particular direction at a particular frequency, we 
mark the point on the appropriate contour where this direction intersects the contour at right 
angles. The gradient indicated by the contour map at this point then gives the group velocity in 
magnitude and direction in accordance with (54). The position vector from the origin to the 
marked point gives the direction of motion of individual wave crests within the pulse and their 
phase propagation constant. This procedure may be carried out analytically by writing the com- 
ponents of (54) parallel and perpendicular to the imposed magnetic field as 
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The derivatives are calculated by differentiating the dispersion equation (55) after converting 
from polar coordinates (/c, P ) to Cartesian coordinates (k\\, k±) in accordance with (19). If we 
reconvert after performing the differentiations, we obtain 
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To calculate the angle g between the direction of motion of a wave packet and the direction of the 
imposed magnetic field we divide (60) by (59), thereby obtaining 
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This equation relates the direction of group propagation to the direction of phase propagation 
and is illustrated graphically in figure 10. Using this relationship we may express the components 
of the group velocity vector in (59) and (60) as functions of g . If the relation between group velocity 
and the angle g is illustrated in polar form, we obtain the diagrams shown in figure 11. In figure 
11 the mirror image in the perpendicular axis should also be drawn and each pattern should be 
rotated round the imposed magnetic field. For a pulse radiated from the origin at time zero each 
diagram in figure 11 gives the radiated pulse pattern at unit time. Pulse patterns at a succession 
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FIGURE 10. Illustrating the relation between the direction of group propagation and the direction 
of phase propagation (a) for the L wave and (b) for the R wave. 
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Figure 11. Illustrating the directional pulse patterns for a source at the origin (a) for the L 

wave and (b) for the R wave. 



of times may be imagined by expanding the scale of the diagrams in figure 11 proportionally to 
elapsed time. 

From figure 11a we see that, for the L wave in its hydromagnetic pass-band, a single pulse 
leaving the source gives rise to two pulses in the radiation field. Moreover, radiation is confined 
to a cone round the imposed magnetic field. The semi-angle of the cone decreases with frequency 
and is given approximately by 
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Thus a source of L waves in the hydromagnetic pass-band constitutes a beam antenna looking 
both ways along the imposed magnetic field, and the beaming becomes more pronounced at 
lower frequencies. As the frequency tends to zero, the pulse pattern shrinks to a pair of points 
moving away from the source in opposite directions along the imposed magnetic field with the 
ionic Alfven velocity. On the other hand, as the frequency tends to the ionic gyrofrequency, the 
group velocity tends to zero. Above the ionic gyrofrequency we are in a stop-band for the L wave 
until a critical frequency is reached at radio frequencies. 

The arrows in figure 11 indicate, in accordance with the angle d p , the direction in which 
individual wave crests move within the pulse. In figure 11a we see that the direction of phase 
propagation varies from longitudinal at the center of the slow pulse to transverse at the tip of 
the fast pulse. At this point group propagation is longitudinal although phase propagation is trans- 
verse. The variation in the magnitude of the phase velocity along the pulse pattern may be seen 
from the dependence of refractive index upon 6 P in figure 7. With the aid of figure 10, this may 
be redrawn to show the variation of refractive index with 6 g as shown in figure 12. At the tip of 
the fast pulse in figure 11a the phase velocity is zero in accordance with the infinity of refractive 
index in figure 12a. 

Figure lib shows the directional pulse patterns for the R wave. At frequencies small com- 
pared with the ionic gyrofrequency this wave is radiated omnidirectionally with the ionic Alfven 
velocity. As the frequency increases, the group velocity along the imposed magnetic field in- 
creases while that perpendicular to the imposed magnetic field decreases. At the bottom of the 
left half of figure lib (/=/*>) the group velocity perpendicular to the imposed field has been reduced 
to zero. At still higher frequencies (/>/<*) radiation is confined to a cone round the imposed mag- 
netic field as shown in the right half of figure lib. Here we again have at the tip of the fast pulse 
a situation in which group propagation is longitudinal while phase propagation is transverse, the 
phase velocity being zero in accordance with the infinity of refractive index in figure 12b. 

Figure 13 shows for group propagation how the hydromagnetic approximation fits onto the 
radio approximation. The left half of figure 13 contains a selection of the directional pulse patterns 
shown in figure lib and is based on the phase propagation diagram in figure 9a. The directional 
pulse patterns in the right half of figure 13 are based on the phase propagation diagram in figure 
9b and therefore on (56). In figure 13 it should be remembered that, in accordance with (53), 
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FIGURE 12. Illustrating the variation of refractive index with direction of group propagation 
(a) for the L wave and (b) for the R wave. 
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FIGURE 13. Illustrating for the R wave the fit between the 
directional pulse patterns derived from (a) the hydro- 
magnetic approximation and (b) the radio approximation. 
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each wave packet moves radially away from the origin. Thus the portions of the pulse that trail 
back to or near the origin are in fact being stretched sideways and are therefore subject to more 
than inverse-square-law divergence. These trailing portions of the pulse are not of great practical 
importance. The trailing portion appearing in the top diagram on the right of figure 13 is A com- 
pletely omitted for the hydromagnetic approximation in the bottom diagram on the left. From 
figure 13 we obtain an integrated picture of how the omnidirectional magnetohydrodynamic wave 
at frequencies small compared with the ionic gyrofrequency becomes transformed into the phe- 
nomena encountered at whistler frequencies. 

The following points should be noted about figure 13: 

(1) Strongest beaming of the R wave along the imposed magnetic field occurs in the radio 
range. As shown by R. L. Smith [1960] the frequency for strongest beaming is approximately 
0.19/v/r and the corresponding half-angle of the cone is about 11°. 

(2) In the frequency range f x <f< x lifMe there is radiated along the imposed magnetic field 
a pulse for which the direction of phase propagation makes a large angle with the direction of the 
imposed field. As shown by Gendrin [1960], the group velocity for this pulse in the frequency 
range fo* < ^f< y^fxie is independent of frequency and is equal to one half the electronic Alfven 
velocity. 

(3) While the transition from omnidirectional propagation to beaming is shown in figure 13 
as occurring near the ionic plasma frequency, in reality it occurs at the frequency / w given by 
(43) for the hydromagnetic approximation and by (45) when an approximate allowance is made 
for the effect of electron inertia at the ionic plasma frequency. 

(4) Because of the effect of divergence in excess of that corresponding to the inverse square 
law, de facto beaming of the radiation occurs below the frequency / x . This is true a fortiori when 
the effect of collisions is taken into account (see sec. 7). Beaming of the R wave along the imposed 
magnetic field is therefore significant for all frequencies large compared with the ionic gyrofre- 
quency and appreciably less than the electronic gyrofrequency. 

It should further be noted that the beaming indicated in figures 11 and 13 arises from the be- 
havior of the magnetoplasma rather than from the characteristics of the source. It would also be 
possible for the source itself to have beaming properties. We may conveniently think of such a 
source as a broadside antenna that emphasizes a particular part of an angular spectrum of plane 
waves. To see which portion of a pulse pattern illustrated in figures 11 and 13 is emphasized by 
beaming of the source, we must look at the directions of phase propagation indicated by the arrows. 
To emphasize a particular portion of a pulse pattern by means of a broadside antenna, the plane of 
the antenna must be at right angles to the corresponding arrow in figures 11 and 13. Thus to 
emphasize the central section of the slow pulse in figure 11a we would orient the plane of the broad- 
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side antenna at right angles to the magnetic axis in the usual way. On the other hand, if we turned 
the plane of the antenna through a right angle, we would then emphasize the tip of the fast pulse; 
energy would slide off the face of the antenna sideways in the direction of the imposed magnetic 
field and in the reverse direction. 



5. Elliptical Polarization Phenomena 

The dispersion relation (21) was obtained by eliminating the quantities (20) from a set of linear 
equations relating them. We can now return to these linear equations and study the ratio of any two 
of the quantities (20). In particular we can study the ratio E y /E x , which describes the elliptical 
polarization of the electric vector in a magnetoionic wave. We have seen that, because of neglect of 
electron inertia in the hydromagnetic approximation, the electric vector vibrates in a plane per- 
pendicular to the imposed magnetic field. Hence the ratio E y /E x gives the complete information 
concerning the elliptical polarization of the electric vector for a magnetionic wave. 

Let us suppose that the linear equations relating the quantities (20) are handled in such a way 
that all these quantities other than the ratio E y /E x are eliminated. We are then left with the 
equations 
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It is the final elimination of E y jE x between (63) and (64) that leads to the dispersion relation (21) 
and hence to (22) for the square of refractive index. If we substitute for the square of refractive 
index from (22) into (63) we obtain for the polarization ratio 
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Alternatively we may eliminate the square of refractive index between (63) and (64) and hence 
obtain a quadratic equation for the polarization ratios of the magnetoionic waves, namely 
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The solution of this quadratic then gives (65). 

The product of the roots of the quadratic (66) has a simple value from which we may derive 
the relation 
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The product of the ratios E y cos 6 P /E X for the two magnetoionic waves is equal to unity. This 
implies a simple relation between the electric ellipses for the two waves that may be described as 
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follows. The two electric ellipses are in a plane perpendicular to the direction of the imposed 
field, and these ellipses may be orthogonally projected onto a plane of constant phase. Either of 
these projected ellipses may be obtained from the other by reversing the direction of rotation and 
turning the ellipse round the direction of phase propagation through an angle V277-. This relation 
between the electric ellipses of the two magnetoionic waves is the same as that familiar at radio 
frequencies [Booker, 1934] . 

The fact that the expression on the right-hand side of (65) is real implies that the components 
of electric field in and perpendicular to the (k, B ) plane vibrate in quadrature. The electric 
ellipse therefore lies in a plane perpendicular to B and has its principal axes in and perpendicular 
to the (k, Bo) plane (see fig. 1). The expression on the right-hand side of (65) is the ratio of the axis 
in the (k, B ) plane to the axis perpendicular to this plane. Positive values imply left-handed 
polarization about B and negative values right-handed polarization. The values ± 1 correspond to 
circular polarization and the values 0, °° to linear polarization. The polarization ratio —jEJE x is 
shown in figure 14 as a function of frequency for a series of values of 6 P . We notice that the polari- 
zation ratios for the L and R waves have opposite signs corresponding to left and right handed 
rotation about B regardless of the direction of propagation. 

From figure 14b we see that the electric ellipse for the R wave at the ionic gyrofrequency is a 
circle perpendicular to the imposed magnetic field regardless of the direction of propagation. 
On the other hand, we see from figure 14a that, for nearly longitudinal propagation of the L wave at 
frequencies small compared with the ionic gyrofrequency, polarization varies rapidly from circular 
to linear. In fact, interchange in the order of the limiting processes p -*0 and /— » effects a 
switch from linear to circular polarization. This is because the curve indicating transition from 
quasi-transverse to quasi-longitudinal propagation passes through the origin in figure 5a. The 



FIGURE 14. Illustrating elliptical polarization properties 
of the electric vector, the magnetic vector and the elec- 
tronic velocity (a) for the L wave and (b) for the R wave. 
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behavior of this curve near the origin is given by 
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from which we see that, when P and co/coa/i are both small, the elliptical polarization of the L wave 
is vitally dependent on the ratio of 0? to co/com*. 

Using figure 10, the angle 6 P between the direction of phase propagation and the direction of 
the imposed magnetic field may be converted into the angle g between the direction of group propa- 
gation and the direction of the imposed magnetic field. In this way figure 14, as well as a number of 
ensuing diagrams, may be replotted in terms of 6 g . While the resulting curves are of interest, 
they are mostly omitted to conserve space. As an alternative, we may remember that 6 P varies 
from to V27T as we move from the center of the slow pulse to the tip of the fast pulse both in figure 
11a and in the right half of figure lib as shown by the arrows. It follows from figure 14 that the 
electric vector is circularly polarized in a plane perpendicular to the magnetic axis at the center 
of the slow pulse, but changes via an elliptical vibration to an almost linear vibration radial to the 
magnetic axis at the tip of the fast pulse. At the tip E vanishes for the hydromagnetic approxi- 
mation although there is some vibration of E along the magnetic axis when electronic inertia is 
taken into account. 

Turning now to the elliptical polarization of the magnetic vector in a magnetoionic wave, 
we remember that, in accordance with (15), this vector vibrates in a plane perpendicular to the 
direction of phase propagation. Nevertheless, the magnetic ellipse is simply related to the electric 
ellipse. From the x and y components of (15) it follows that 

#£ = _£j> . (69) 

Hy E X 

This equation implies that the orthogonal projection of the magnetic ellipse onto the plane of the 
electric ellipse is identical in shape with the electric ellipse except for rotation through Sin round 
the direction of the imposed magnetic field. Figure 14 may also be interpreted therefore to give 
the elliptical polarization of the magnetic field of a magnetoionic wave. 

The electronic ellipses of the magnetoionic waves may be studied by using the linear equations 
relating the quantities (20) to obtain the ratios of v ex , v ey and v ez - It turns out that the plane of the 
electronic ellipse does not in general coincide with either the plane of the electric ellipse or the 
plane of the magnetic ellipse. However the projection of the electronic ellipse onto a plane per- 
pendicular to the imposed magnetic field is identical in shape and orientation to the projection of 
the magnetic ellipse. Hence the left hand side of (65) may be replaced by 

- ^y .£*£ • Vex 

E x H y V ey *■ ' 

and figure 14 may be interpreted accordingly. However the normal to the electronic ellipse is 
turned away from the vector B towards the vector k through the angle 

0e = -tan- 1 — (71) 

Vey 

and it follows from the linear equations relating the quantitites (20) that 

Vez a . n COtoMi (kc\* ( . E y \ 

— = cos0psin0p — £"( — ~ir ' (72) 

Vey P (0%i \Q)J \ E X J ^ 

Using (71) and (72) we may convert figures 7 and 14 into figure 15. 
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FIGURE 15. Illustrating variation of the orientation of the 
plane of the electronic ellipse (a) for the L wave and (b) 
for the R wave. 



FIGURE 16. Illustrating elliptical polarization properties 
of the electric current vector and the ionic velocity (a) 
for the L wave and (b) for the R wave. 



For the L wave in its hydromagnetic pass-band figure 15a shows that the values of e are nega- 
tive. This implies that, as the direction of phase propagation turns away from the direction of the 
imposed magnetic field, the normal to the electronic ellipse turns in the opposite direction. How- 
ever, these are circumstances in which great importance attaches to the motion of the ions even 
from the viewpoint of electric current. 

For the R wave we have positive ordinates in figure 15b, showing that the normal to the elec- 
tronic ellipse turns away from B in the same direction as k. Moreover, for f>f x , the electronic 
ellipse does not depart radically from perpendicularity to the direction of phase propagation. In 
the right half of figure lib (/*>/*>) the electronic motion is circularly polarized in a plane perpen- 
dicular to the magnetic axis at the center of the slow pulse, but changes via an elliptical motion to 
an almost rotary vibration round the magnetic axis near the tip of the fast pulse. At the tip the 
electronic motion vanishes for the hydromagnetic approximation, although there is some motion 
along the magnetic axis when electronic inertia is taken into account. On the other hand, for 
/</«„ figure 15b shows that, for the R wave, 6 e at first increases from zero with 6 P and then de- 
creases again to zero as 9 P increases further to ¥2 it. Thus, in the left half of figure lib (/</»), 
the plane of the electronic ellipse is perpendicular to the magnetic axis where the pulse pattern 
crosses this axis, begins to turn the same way as the direction of phase propagation as we move 
round the pulse pattern, but returns to a position normal to the imposed magnetic field where the 
pulse pattern crosses the perpendicular axis. For frequencies small compared with the ionic gyro- 
frequency, the plane of the electronic ellipse for the R wave does not depart much from perpendi- 
cularity to the magnetic axis; the electronic motion varies from circular polarization where the pulse 
pattern crosses the magnetic axis to almost linear vibration at right angles to B in the (k, B ) 
plane where the pulse pattern crosses the perpendicular axis. 
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For (65) the quasi-longitudinal and quasi-transverse approximations are 

_ -E_y_ •#£_ -Ve£ , 1 

J E X J H y - J v ey ~-cos6 p (QL) (73) 

E x My V ey (O \ 0J 2 vj 0)%^ P ^ K'*) 

CO 

_.Ey__.H J r_.Vex__ (OM 1 /^yp y^ 

7 £, J H y ~ J v ey 1+ o^_^sin 2 p l ^ ' Aj * 

w m ^m (75) 

A number of features of figure 14 may be traced with the aid of these approximations. 

For the hydromagnetic approximation, the ionic ellipses for the magnetoionic waves lie in 
planes perpendicular to the imposed magnetic field because, as already mentioned in section 2, 
there is no vibration of ions parallel to the imposed magnetic field, there being no electric field 
in this direction to drive them. On the other hand, the ionic ellipses are not in general identical 
in shape with any of the ellipses so far discussed. From the linear equations relating the quantities 
(20) it follows that 

— +(-jZ* 

.Vjx_ (OMi \ E x , 

J v lv 1+J o_(_ 4 E JL \' < 76 > 



a>Mi \ Ej 

With the aid of this equation curves for the polarization ratio for the ions may be deduced from 
figure 14 giving the results shown in figure 16. 

From figure 16a we see that, for the L wave at the ionic gyrofrequency, the ions are circularly 
polarized in a plane perpendicular to the imposed magnetic field regardless of the direction of 
propagation. We are here at the high frequency end of the hydromagnetic pass-band for the L 
wave where the wave is known as the ion cyclotron wave. At the center of the slow pulse in figure 
11a the ionic polarization is circular. On the other hand, at the tip of the fast pulse, the ionic 
motion involves an oscillation around the magnetic axis combined with a quadrature oscillation 
radial to the magnetic axis. 

From figure 16b we see that, for the R wave, the sign of the ionic polarization ratio for trans- 
verse propagation is not the same as for longitudinal propagation. This means that, when the 
direction of phase propagation in the R wave is sufficiently transverse to the imposed magnetic 
field, the ions execute ellipses in planes perpendicular to the imposed magnetic field in the left- 
hand sense. Thus, whereas for the hydromagnetic approximation the magnetoionic waves may 
be classified into left-handed and right-handed waves from the viewpoint of electromagnetic theory, 
the same is not true from the viewpoint of fluid mechanics. However, when the direction of 
rotation of the ions for the R wave is reversed, they are normally executing quite narrow ellipses, 
and the ionic motion is almost a linear vibration perpendicular to B in the (k, B ) plane. In the 
left half of figure lib (/</ x ) the fluid motion varies from circularly polarized near the magnetic 
axis to almost linear vibration in the direction of propagation near the perpendicular axis. At 
the tip of the fast pulse on the right half of figure lib (f>f x ) the fluid vibration is mainly radial 
to the magnetic axis, although there is some quadrature rotary vibration round the magnetic axis. 

The electric current density J is derivable from the ionic and electronic vibrations in accord- 
ance with (6). Since the plane of the electronic ellipse is tilted with respect to the plane of the 
ionic ellipse, the same is true for the current ellipse. However, it emerges that 



Jy = Vix 
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Figure 17. Illustrating variation of the orientation of the 
plane of the current ellipse for the R wave. 
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so that the orthogonal projection of the current ellipse onto a plane perpendicular to the imposed 
magnetic field is identical with the ionic ellipse except for rotation through a right angle. Thus 
figure 16 may also be interpreted as giving the projection of the current ellipse onto a plane per- 
pendicular to the imposed magnetic field. However, the normal to the current ellipse in general 
makes with the direction of the imposed magnetic field a nonzero angle 6j measured from B () 
towards k and given by 



0,/ = — tan l — • 



(78) 



From (6), together with the linear equations relating the quantities (20), we may deduce that 



tan 0./ = tan 



1 + 






J E x 



CO Mi \ Ex) 



(79) 



Using this equation we may calculate 0j from the information plotted in figure 14. However, for 
the L wave in its hydromagnetic pass-band, the quantity in square brackets in (79) is practically 
unity if com? ^ o*m, so that the current ellipse is almost exactly perpendicular to the direction of 
phase propagation. For the R wave, however, departures from this occur as shown in figure 17. 
These departures are particularly noticeable for/</ x and 6 P > tan -1 fmlfMi- For f<f x the angle 
6 j increases from to about V2 77 as 6 P increases from to about tan -1 fmlfiau but then increases 
further to it as 6 P increases to V2 it. During the latter stage, therefore, the current vector rotates 
left-handed about the direction of the imposed magnetic field even though we are discussing the R 
wave. This is a consequence of the corresponding behavior of the ionic ellipse. Thus, in the left 
half of figure lib (/</oo), the current ellipse is perpendicular to the direction of the imposed mag- 
netic field both where the pulse pattern crosses the magnetic axis and where it crosses the per- 
pendicular axis, and the plane of the ellipse makes a complete half revolution between the two 
positions. Moreover, we see from figure 17 that this is true even in the limit of zero frequency. 
A number of features of figure 16 may be seen by means of the quasi-longitudinal and the quasi- 
transverse approximations obtained by substituting from (73), (74), and (75) into (76). We obtain 





CO + 1 


. Jy . Vix 


CO\fi COS dp 


Jx Viy 


l+ w \ 



(80) 
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, -=-+ 1+- i r tan 2 P 

y - ^ " 1+ 1+^-^ Van 2 ., 

, 1 — 11 + — — sin 2 6 P 

■ J z- J « y -^ ri& *\ — (QT ' X) - (82) 



For strictly transverse propagation, (81) and (82) become 



• Jy . Vix (*>(t>Mi 1 m 

- j z =J ^ = ^ — ^" (T ' X) (84) 



and from these equations the behavior of the 90° curve in figure 16b may be verified, thus con- 
firming in these circumstances the left-handed behavior of the R wave so far as electric current 
and fluid motion are concerned. 



6. Comparison of Ionic and Electronic Motions 

To compare the ionic and electronic motions in a magnetoionic wave it is appropriate not 
merely to compare their elliptical polarizations, but also to compare the actual amplitudes of 
vibration. From the linear equations relating the quantities (20), we may derive expressions for 
the ratio of electronic velocity to ionic velocity (a) for the component of motion perpendicular 
to the (k, Bo) plane and (b) for the component of motion in the (k, B ) plane perpendicular to B . 
We obtain 

Vex _ -> 0)/(OMi 



Vix JVixlViy (85) 

and 

V ey -. CO ( .Vix\ 

= 1 1/ - )' (86) 

Viy ixiMi \ ViyJ V "^ 

For the component along the imposed magnetic field the ionic motion is zero for the hydromagnetic 
approximation. By means of (85) and (86) we may convert figure 16a into figure 18. 

We notice that, for longitudinal propagation of the L wave at frequencies small compared 
with the ionic gyrofrequency, the ionic and electronic motions are almost indetical so that the 
corresponding contributions to the electric current almost cancel. It should be noted, however, 
that it is the existence of the "small" difference current that causes the change in the velocity 
of propagation from the free space value to the Alfven value, and this change is normally far from 
"small." 

For the L wave at the ionic gyrofrequency, there is no motion of the electrons in a plane per- 
pendicular to the imposed magnetic field regardless of the direction of propagation. However, 
it follows from figure 15a that there is motion of the electrons along B except when the direction 
of phase propagation is also along B . Thus, for the ion cyclotron wave, there is circular motion 
of the ions in planes perpendicular to the imposed magnetic field combined with linear vibration 

488 




0.2 0.4 0.6 0.8 1.0 1.2 0.2 0.4 0.6 0.8 1.0 
RATIO OF FREQUENCY TO IONIC GYRO FREQUENCY 



FIGURE 18. Illustrating for the L wave variation of the 
ratio of electronic velocity to ionic velocity (a) for the 
component of motion perpendicular to the (k, Bo) plane 
and (b) for the component of motion in the (k, B ) 
plane perpendicular to B . 
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Figure 19. Illustrating for the R wave variation of the 
ratio of electronic velocity to ionic velocity (a) for the 
component of motion perpendicular to the (k, B ) plane 
and (b) for the component of motion in the (k, B ) 
plane perpendicular to B . 



of electrons along the imposed magnetic field, the latter motion vanishing for strictly longitudinal 
propagation. 

In figure 11a the ionic and electronic motions are circularly polarized and in phase at the center 
of the slow pulse, with the positive contribution of the ions to the electric current outbalancing 
the negative contribution of the electrons. As we move to the tip of the fast pulse, the ionic ellipse 
remains perpendicular to B whereas the electronic ellipse rotates in the opposite direction to k 
as shown in figure 15a. The combination produces a current ellipse that turns so as to be almost 
exactly perpendicular to the direction of phase propagation. At the tip of the fast pulse the elec- 
tric current is dominated by oscillation of the electrons along the imposed magnetic field. 

For the R wave we may use (85) and (86) to convert figure 16b into figure 19. We notice that, 
at frequencies large compared with the ionic gyrofrequency, the electronic motion far exceeds 
the ionic motion. This is what, at sufficient high frequencies, justifies use of the radio approxi- 
mation. However, it should be noted from figure 19b that v e y < v w when the direction of phase 
propagation is nearly transverse to the imposed magnetic field. In these circumstances, there- 
fore, the negative contribution to the electric current from the small component of electronic 
motion in the (k, B ) plane perpendicular to B is outbalanced by the positive contribution from 
the ionic motion. This is why, for nearly transverse propagation, the electric current has left- 
handed polarization for the R wave even though the electronic motion has right-handed polariza- 
tion (compare fig. 16b with fig. 14b). 

Equations (12) and (13) show that, except for propagation strictly along the imposed magnetic 
field, the magnetoionic waves involve compressions and rarifications in both electron density 
and ion density. In fact these equations state that the fractional increase in density is equal to 
the product of the refractive index and the component of particle velocity in the direction of phase 
propagation, the velocity being expressed as a ratio to the standard velocity used in specifying 
refractive index (either the free space velocity or the ionic Alfven velocity). Consequently, if 
the ratio to the ionic Alfven velocity of the particle velocity in the direction of phase propagation 
is known, we merely multiply this by the appropriate ordinate in figure 7 to obtain the fractional 
increment in particle density. Near the infinite ordinates in figure 7 the assumption of a linear- 
ized cold plasma breaks down. 

It is also interesting to compare the fluctuations in electron density with those in ion density. 
From (12) and (13) we obtain 



8Ne _ Vey sin Op + Vez cos 6p 

8Nj Vjy sin Op 
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Figure 20. Illustrating the variation in the ratio of 
electron density to ion density for the R wave. 
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and using (71), this may be rewritten 



We 
Wi' 



tan 6 e 
tan Or 



(88) 



We may therefore use the information plotted in figures 15, 18b, and 19b to derive 3N e l8Ni. For 
the L wave in its hydromagnetic pass-band this ratio is almost exactly unity independent of the 
direction of propagation. Except for strictly longitudinal propagation, there are fluctuations in 
both the ion density and the electron density but very little fluctuation in charge density. Except for 
frequencies small compared with the ionic gyrofrequency, the motions of the ions and electrons are 
different, but the components of motion perpendicular to the direction of phase propagation are al- 
most exactly the same so that the net fluctuation in space-charge is practically zero. For the R wave, 
on the other hand, 8N e /8Ni departs from unity for frequencies large compared with the ionic 
gyrofrequency. The dependence on direction of propagation is relatively slight as shown in figure 
20, but it should be remembered that both 8N e and 57V; vanish for strictly longitudinal propagation. 
Motion of the ions plays an essential role in the fluctuation of space-charge until the frequency 
is high compared with/oo. 



7. Ionic Attenuation 

It is well known that, for the radio approximation, it is easy to allow for electronic collisions 
by replacing the electronic mass m e by m e U e where 



U e =l—J 

(x) 



(89) 



In the same way, for the hydromagnetic approximation, it is easy to allow for collisions of ions 
with neutral particles provided that the neutral gas is sufficiently massive to take no significant 
part in the wave motion. If v\ is the ionic collisional frequency involved, we add the term —miVi\i to 
the right hand side of (2). In (11) this has the effect of replacing the ionic mass m\ by rmUi where 






(90) 



Consequently, to allow for ionic collisions with a neutral gas that takes no significant part in the 
wave motion, we modify the dispersion relation (21), and all subsequent equations, by replacing 



(o M i by (OMilUi 
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FIGURE 21. Illustrating the variation in the rate of ionic 
attenuation (a) for the L wave and (b) for the R wave. 
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and 



<*>jvi by (a 2 Ni IUi. 



(92) 



To calculate the first order effect of ionic attentuation we may make the replacements (91) 
and (92) in the dispersion relation (55) and then expand to the first order in powers of (U% — 1). We 
obtain in this way 



1 co \dUi/Ui=i 



(93) 



so that the rate of attenuation in nepers per unit distance is 

= Pj (dk\ 
co XdUju^i . 



(94) 



By differentiating the modified dispersion equation with respect to U\ and then putting U, = 1, we 
derive 
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(95) 
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This is the number of nepers of ionic attenuation suffered in the distance Aijvu which is the wave- 
length obtained by converting the ionic collisional frequency into wavelength using the ionic 
Alfven velocity. The rate of attenuation derived from (95) is illustrated in figure 21. From figure 
21a we see that the attenuation of the L wave becomes large as the frequency increases to the 
ionic gyrofrequency, where the wave is the ion cyclotron wave. From figure 21b we see that, 
for the R wave, the attenuation at low frequencies is half a neper per distance A\\v[ and decreases 
with increasing frequency except when the direction of phase propagation is nearly transverse 
to the imposed magnetic field. For the latter situation there is a maximum of attenuation near 
the frequency f x given by (43) for the hydromagnetic approximation and by (45) when electron 
inertia is taken approximately into account at the ionic plasma frequency. If attenuation is im- 
portant, it is more important for propagation perpendicular to the magnetic axis than for propaga- 
tion along it. Thus the effect of attenuation in the left half of figure lib is further to suppress 
those features of the pulse pattern that are subject to more than inverse-square-law divergence. 
Ionic attenuation also has some effect on the elliptical polarization properties of the magneto- 
ionic waves. In particular the interpretation of (67) is affected. The projections of the electric 
ellipses for the two magnetoionic waves on a plane perpendicular to the direction of phase prop- 
agation, instead of being perpendicularly located, are turned through equal angles in opposite 
directions. Also, the condition for equality of the refractive indices of the two magnetoionic 
waves, instead of occurring for the hydromagnetic approximation when the direction of phase 
propagation is strictly transverse to the magnetic field, occurs at a small angle to this direction 
equal to Vi(x)Mil(o 2 Ni . 
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